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Abstract 

Recent models of soft diffraction include a hard pomeron pole besides the 
usual soft term. Such models violate the black-disk limit around Tevatron 
energies, so that they need to be supplemented by a unitarisation scheme. 
Several such schemes are considered in this letter, where we show that 
they lead to a large uncertainty at the LHC. We also examine the impact 
of unitarisation on various small- i observables, the slope in t of the elastic 
cross section, or the ratio of the real to imaginary parts of the scattering 
amplitude, leading to the conclusion that the existence of a hard pomeron 
in soft scattering may be confirmed by LHC data. 

PACS numbers: 13.85.Lg, 13.85.-t, 13.85.Dz, 12.40.Nn, 11.80.Fv 
Keywords: Saturation, unitarisation, hard pomeron 



Introduction: hard poles 

Experimental data reveal that total cross sections grow with energy. This means 
that the leading contribution in the high-energy limit is given by the rightmost 
singularities in the complex- j plane, the pomerons, with intercepts exceeding 
unity. In the framework of perturbative QCD, the leading singularity is expected 
to exceed unity by an amount proportional to a s pp. At leading-log, one obtains 
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parameter 


value 


error 


parameter 


value 


error 


a s (0) 


1.0728 


0.0008 


a H (0) 


1.45 


0.01 


Sp 


56.2 


0.3 


Hp 


0.10 


0.02 




32.7 


0.2 


H n 


0.28 


0.03 


Sk 


28.3 


0.2 


H K 


0.30 


0.03 




0.174 


0.002 


Hj 


0.0006 


0.0002 



Table 1: Parameters of the leading singularities of the fits of ref. [2] for y/s from 
5 to 100 GeV. 

a leading singularity at J — 1 = 12a s log2/V. From analytic S- matrix theory, it 
is expected that the singularities of the exchanges do not depend on the external 
kinematics. Hence, although perturbative calculations can be justified only for 
the scattering of far off-shell particles, the singularities that they predict should 
remain valid in the long-distance range. 

In a recent study j2], we have indeed found that forward data (total cross 
sections and the ratios of the real part to the imaginary part of the amplitude) 
could be fitted well by a combination of a soft pomeron (which would be purely 
non perturbative) and a hard pomeron. The expression of the leading terms of 
the total cross sections for the scattering of a on p becomes 

eg, = ^f raA (^) (1) 

with P the beam momentum in the target frame, s, t and u the Mandelstam 
variables, and A the hadronic amplitude 

ZmA(s) « H a (-J + S a (-J (2) 

with S\ = 1 GeV 2 , and the parameters given in Table 1. The inclusion of these 
two pomerons, together with the use of integral dispersion relations, and the 
addition of sub-leading meson trajectories {p/uj and a/ f), leads to a successful 
description of all pp, pp, ^p, K ± p, jp and 77 data for ^/s < 100 GeV. 

The problem however is that such a fit does not extend to high energies. 
Indeed, the fast growing hard pomeron leads to a violation of unitarity (an elastic 
cross section bigger than the total cross section) for values of s smaller than 1 
TeV. 

In [2], we used a simple ansatz to unitarise the hard pole. However, we did 
not examine the uncertainties linked to it, as we simply wanted to show that it 
was possible to extend the fit to higher energies. In this letter, we propose to 
reconsider this question and to give an estimate of the uncertainties on total cross 
section measurements at the LHC. 
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Section 1 will be devoted to a reminder of unitarisation, and of the impact- 
parameter (b) formalism. Section 2 will consider the "minimal" unitarisation, 
which cuts off the amplitude in impact parameter once it reaches the black-disk 
limit. Section 3 will consider analytic unitarisation schemes. Finally, putting 
everything together, we shall show that the cross section at the LHC should be 
large, of the order of 150 mb. 



1 Unitarity and black-disk limit 

The simplest expression of the unitarity constraints can be obtained from partial- 
wave amplitudes. We can write 

A(M) = 87r£(2Z + l)^(s)P;(cos0 s ). (3) 
i 

At high energy and small angle, one can rewrite I ~ so that the partial- wave 
decomposition can be rewritten in impact-parameter space as 

A(s, t = -q 2 ) = 2sJ d 2 be- i& T{s, b) (4) 

In terms of J-{s, b) = J d 2 qe z ^ b A(s,t)/ (8n 2 s), the total and elastic cross sections 
are given by 

cr to t = An J bdb^smF^s, b) (5) 

and 

a el = 2n Jbdb\F(s,b)\ 2 . (6) 

One can then show that unitarity of the S matrix, SS^ = 1, together with 
analyticity, and with the normalisation used in (JS]) and (jH]), requires that 

< \F(s, b) | 2 < 2%mF(s, b) < 4. (7) 

At high energies, as SmjF oc s A with A > 0, the scattering amplitude F(s, b) 
reaches the unitarity bound for some value b u (s) of the impact parameter. 

However, before this happens, another regime is reached: one gets the maxi- 
mum inelasticity if 23mjF(s, b) — \ J r (s,b)\ 2 is maximum, i.e. if 

SmF(s,b) = l. (8) 

This is the region where the proton becomes black, and it is usually referred to 
as the black-disk limit. The imaginary part of T is usually noted as T(s, b) and 
called the profile function. Saturation usually refers to the black-disk limit. We 
shall use it to mean that we reach the maximum possible inelastic amplitude at 
some distance bs (s) between scattering particles. 
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In QCD, this can be thought of as a consequence of the growth of the gluon 
density at small x, which must be tamed by non-linear effects connected to the 
next-to-leading terms. One needs to note that the role of non-perturbative ef- 
fects connected with confinement [3] and with large impact parameters makes 
the analogy with perturbative saturation questionable, as in the present work 
saturation will first occur at small b. 

One expects saturation to tame the growth of a tat . For example, in the loop- 
loop correlation model (LLCM) [4] based on the functional integral approach to 
high-energy collisions [5] , the T- matrix element for elastic proton-proton scatter- 
ing reads 

A pp (s,t) oc is J d 2 be icl *T(s,b) 

T(s,b) = Jdz 1 d 2 r 1 Jdz 2 d 2 r 2 \i>p(z 1 ,fi)\ 2 \'ilj p (z2 1 r2)\ 2 

x 1 - S DD (s,b, z 1 ,r 1 ,z 2 ,r 2 ) (9) 

with ip the proton light-cone wave function and Sdd the dipole-dipole S matrix 
element. The saturation of the profile function is a direct consequence of S- 
matrix unitarity. In this model, the saturation regime at small b is reached only 
at very small x ~ 10~ 10 and very high energies y/s > 10 6 GeV, where there is 
a transition from a power-like to an log 2 -increase of cr^(s), which then respects 
the Froissart bound [6]. A similar result, in the framework of the dipole picture 
of soft processes, also leads to a taming a tot but at lower energy [7]. 

If we take a single simple pole for the scattering amplitude, with an exponen- 
tial form factor of slope d, 

A(s,t) = S p s as{0)+a,t e dt (10) 

then the radius of saturation and its dependence on energy can be obtained 
analytically: 

6 s(s )^ W log,)log( ^^ ) (11) 

Approximating a tot w 27rfo|(s), one sees that the total cross section grows loga- 
rithmically at medium energies and like log 2 s at very high energies. 

In our model [21 [8], the pp elastic scattering amplitude is proportional to the 
hadrons form-factors and can be approximated at small t as: 

A{s,t) = [H p F H (t)(s/ Sl ) a « (0) e< * log(s/si) 

+ S p F s (t) (s/si) Qs(0) e a 's * log(s/si) ] (12) 

where the couplings and intercepts are given in Table 1. The study of ref. [S] of 
small-t elastic scattering shows that the slope a' s of the soft pomeron trajectory 
is slightly higher than its classical value [9j[10], and we shall take a' s = 0.3 GeV 2 . 
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parameter 


value 


parameter 


value (GeV~ 2 ) 


h 


0.55 


d\ 


5.5 


h 2 


0.25 


d 2 


4.1 


h 


0.20 


d 3 


1.2 



Table 2: Parameters of the elastic pomeron form factor, see Eq. (fT3j) . 

The slope of the hard pomeron trajectory is evaluated [8j[TT| to be a! H = 0.1 GeV 2 . 
The normalisation Si = 1 GeV 2 will be dropped below and s also contains im- 
plicitly the phase factor exp(— in/2), corresponding to crossing symmetry. 

A small-t analysis [8] indicates that the form factor F$(t) is close to the square 
of the Dirac elastic form factor, and can be approximated by the sum of three 
exponentials [ID] . 

= ( im l-2.m i y 



V 4m 2 -t l-t/A 2 ) 
« h x e dl 1 + /i 2 e d2 ' + h 3 e d3 *. (13) 

where m p is the mass of the proton, A 2 = 0.71 GeV 2 . The other parameters are 
given in Table 2. For the hard pomeron, the form factor is rather uncertain [8], 
and we assume that it can be taken equal to Fs(t). 

We then obtain in the impact parameter representation a specific form for the 
amplitude in b space, !Fq(s, b) [12], which we show in Fig. CQ 

^oM) = s a ^ exp(-b 2 /rls) + (S ^ H) 

with r 2 5 = 4 {di + a' s log(s)) (14) 
r 2 H = 4(4 + 04 log(*))- ( 15 ) 

One can see that at some energy and at small b, T (s,b) = QmFo(s,&) reaches 
the black disk limit. For our model, this will be in the region y/s ~ 1.5 TeV. 

What happens then is largely unclear. There is no unique procedure to sat- 
urate the amplitude, and saturation and unitarisation will be important issues 
at the LHC. In the following, we shall consider two main schemes: the minimal 
one, where the saturation freezes the profile function at 1: this is the minimum 
possibility for the restoration of unitarity, as it does not affect in any way the 
low-energy data. In the next section, we shall also consider the predictions of an 
eikonal scheme. 

Saturation of the profile function will surely control the behaviour of o t ot 
at higher energies. We assume that once it reaches 1, the amplitude does not 
change anymore and remains equal to i: recombination must be maximal for 
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Figure 1: The profile function for proton-proton scattering: the hard line and 
circles - at \fs = 2 TeV without and with the saturation; the dashed line and 
squares - at y/s = 14 TeV without and with saturation, the dash-dotted line - 
the eikonal form (pDI) at */s = 14 TeV. 

black protons. But this freezing of the profile function must be implemented 
carefully: one cannot simply cut the profile function sharply as this would lead 
to a non-analytic amplitude, and to specific diflractive patterns in the total cross 
section and in the slope of the differential cross sections. Furthermore, we have 
to match at large impact parameter the behaviour of the unsaturated profile 
function. 

We use an analytic interpolating function, which is equal to one for large 
impact parameters and which forces the profile function to approach 1 at the 
saturation scale b s as a Gaussian. Analyticity of the function enables us to use a 
complex s as before to obtain the real part. 

We assume saturation starts at a point b a little before b s , and that the profile 
function is 1 for b < b . The saturated profile, T s (s, b), is otherwise given by 



We find that we need to assume that the scale in the Gaussian is s-dependent 
because the slope of the profile function decreases with energy. A reasonable 
match is provided by b y = 32/ logs, which is about 5 GeV -1 at the Tevatron. 
We shall give our results for b 2 z = 2 GeV~ 2 and b /b s = 97.5%. 



r s (s,b) 




(16) 



l + (r( S ,0)-l)exp[-((^)) 2 
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Figure 2: The total cross section as a function of yfs, for the bare amplitude 
(short dashes), the saturated amplitude (plain curve), the eikonalised amplitude 
(dash- dot- dot), and for a renormalised eikonal (long dashes). 

We show in Fig. [2] the behaviour of the total cross section at high energies. 
We see that saturation brings in a significant decrease of the LHC cross section. 
However, it is also clear that the simple saturation considered here is not enough, 
as the total cross section at the Tevatron will be 85 mb, which is 2 standard devi- 
ations from the CDF result. It may be that our previous estimate of the pomeron 
coupling or intercept is too big. We show in Fig. [3] the effect of changing either. 
We see that if saturation is the driving mechanism, then the coupling would need 
to be reduced^ to 0.06 or the intercept to 1.41 in order to accommodate the Teva- 
tron point in this scheme. As we want to show qualitatively what the effect of 
saturation might be, we keep the parameters of Table 1 in the following. 

The saturation regime will have some major effects at the LHC. We show in 
Fig.H]that the elastic cross section will be somewhat affected, and that its growth 
will be tamed: the ratio <y e i/u tot will start a slow growth towards 0.5. But more 
importantly, the small- 1 data will look quite different. We show in Fig. [5] the 
behaviour of the ratio of the real-to-imaginary part of the cross section, both in 
the bare and in the saturated case. From it, we see that the small-t slope of p 
will be one of the most striking features of saturation. If this is not measurable 
at the LHC, then one can also consider the slope of the differential elastic cross 
section, which we show in Fig. [6j We see that saturation increases the slope at 
small t, and predicts a fast drop around \t\ = 0.25 GeV 2 , when one enters the 

1 Note that this is still four times larger than the value advocated in [T3]. 
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Figure 3: The total cross section at the Tevatron (lower curves) and at the LHC 
(upper curves) as a function of the intercept (left) and of the coupling (right) of 
the hard pomeron for a saturated amplitude (plain curves and short dashes) and 
for an eikonalised amplitude (long dashes and dash-dots) . 




102 103 104 102 103 104 
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Figure 4: The elastic cross section as a function of s (left), and its ratio to the 
total cross section (right), for the bare amplitude (short dashes), the saturated 
amplitude (plain curve), the eikonalised amplitude (dash-dot-dot), and for an 
extended eikonal (long dashes). 
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Figure 5: The ratio of the real to the imaginary part of the amplitude as a function 
of t, for the bare and the saturated amplitudes at various energies: 100 GeV (plain 
curve), 500 GeV (long dashes), 5 TeV (short dashes) and 14 TeV (dash-dotted 
curve). 

region of the dip. 

In fact, saturation naturally predicts a small increase of the slope with t at 
small t. To understand this, let us take the simple form of the black disk with a 
sharp edge at radius R. The scattering amplitude can then be represented as 

In this case, the slope of the differential cross section at small momentum transfer 
will be 

B BDL ~R 2 /A +#732 |t|. (17) 

Hence the slope will grow with increasing \t\ at small momentum transfer, as can 
be seen in Fig. El 



2 Analytic unitarisation schemes 

As we have seen, in b space, unitarisation can be written 

2%mF - (SmJ^) 2 - (KeJ* 7 ) 2 = g in > 0. (18) 

The general solution can be obtained iff g>j n < 1, and can be written T = i(l — 
(1 — (7m)e 1 *). We can rewrite it using the opacity Q so that (1 — g in ) = e~ n 

T = e- n+i *) =i(l- e^ s ' b) ) . (19) 
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Figure 6: The slope of the elastic differential cross section as a function of t, for 
the bare and saturated amplitudes at various energies: 100 GeV (plain curve), 
500 GeV (long dashes), 5 TeV (short dashes) and 14 TeV (dash-dotted curve). 

Any unitarisation method has to lead to such a form for the amplitude. The 
ambiguity however comes when one tries to identify gi n and $ in formula (fT9l) 
with the physics input. It is known from potential models that in non-relativistic 
physics one can think of the Taylor expansion of ffl9|) as a description of successive 
interactions with the potential. Here, however, we have no potential, so that the 
identification of each term with successive pomeron exchanges is not obvious. 

The usual approach is to assume that the Taylor expansion of (|T9j) is such 
that the n th term corresponds to n-pomeron exchange. In this case, we take the 
eikonal form for the scattering amplitude. 

A e (s,t = 0) = 2 /d 2 6 [1 - exp(iF (>, 6))] . 

Before giving the results in this approach, a few comments are in order. First 
of all, the eikonal is only a model. Indeed, it is known [TJj that it does not 
reproduce properly the s-channel cuts of the scattering amplitude coming from 
multiple exchanges. For instance, already the second term, corresponding to the 
two-pomeron cut, could have a suppression due to the structure of the proton 

Furthermore, the eikonal does not always guarantee unitarisation. Although 
it obeys Eq. ([7]) for all values of b, it can produce, after integration, amplitudes 
that violate the Froissart bound [5J, depending on the dependence of the form 
factor. To see this, take the eikonal in factorised form 

X (s,b) = h(s) f(b) } (20) 

with h(s) = s A , and assume simple functional forms for the form factor f(b), 
which allow an analytical treatment. If one considers a Gaussian form 

f(b) ~ exp(-6 2 / J R 2 ), (21) 
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one obtains 

A(s,t = 0) ~ i R 2 (r(0,s A ) + 7 + A logs), (22) 

where 

roc 

r(a,z) = / t^ 1 e~ l dt 



and, in our case, 

T(0,s A ) -»• 0, s -> oo. (23) 

If i? 2 is independent from s, we have 

<r toi ~ log(s). (24) 

whereas for R 2 growing like log(s), we obtain 

a tot ~ log 2 (s), (25) 

so the total cross section respects the Froissart bound. However, if we take a 
polynomial form, such as that resulting from dipole-dipole interactions [15j, 



we obtain 



f(b) ~ \, (26) 



A(s,t)~i / [1 - eX p(-s A 2 /)] = 2v^Fs A/2 . (27) 

Jo V\fv 

with y = l/b A . So, in this case, the scattering amplitude does not satisfy the 
Froissart unitarity bound. This violation is not due to the divergence of f(b) 
at small b. Indeed, if we introduce an additional small constant radius r which 
removes the singular point b = in f(b) and take 

m - ^ (28) 

the answer, after some complicated algebra, is 

A(s,t = 0) ~ i L[ 7r ^ exp [-^/(2r 4 )] [J ( S A /(2r 4 )) + / 1 ( S A /(2r 4 ))]. (29) 
The asymptotic value of the Modified Bessel functions is 

J ,i(s A /(2r 4 )) ~ . (30) 



x /vrs A / 2 ' 
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Hence we again obtain for asymptotic high energies 



AM)~i^=— (31) 

which again does not obey the Froissart bound. 

So, in the following, we shall use the eikonal as a simple example. The form 
factors that we use lead to a cross section that does respect the Froissart bound. 
Furthermore, we shall be able to see whether the saturation effects that we found 
depend on the picture of saturation, or depend only on the onset of unitarising 
cuts. 

We show in Fig. [I] that, in this case, the profile function never saturates: the 
cuts actually reduce the cross section from the start, and never allow F(s, b) to 
become larger than 1. This means, as is shown in Fig. [2] that eikonalisation will 
give a suppression, even at lower energies. So one needs to modify slightly the 
parameters of Table 1, to recover the low-energy fit. We find that multiplying 
the couplings by 1.2 provides such an agreement, and refer to this as the "renor- 
malised eikonal" , shown in Fig. [2j We see that both curves are very close at the 
LHC, but exceed considerably previous estimates [16]: whereas the total cross 

+4.1 

section at the LHC was predicted to be 111.5 ± 1.2 mb, it is now 152 mb. 

In a similar manner, Fig. H] shows that the elastic cross section also gets large 
corrections from the cuts at lower energies. Again, the renormalised eikonal 
is close to the saturated curve at low energy, but this time deviates at higher 
energies, so that the ratio of the elastic cross section to the total cross section 
changes by about 10%. 

We had seen that a striking feature of saturation was the behaviour of the real 
part of the amplitude at small t, as the p parameter would have a drastic change 
of slope. We find that such an effect will also be present in the eikonal case (see 
Fig. El), although the real part will be much bigger in this case. Similarly, we also 
show in Fig. [7] the behaviour of the slope of the differential elastic cross section. 
Again, we see that it increases at small t, and then decreases towards the dip. 



3 Conclusion 

Now the structure of the diffractive scattering amplitude cannot be obtained 
from first principles or from QCD. The procedure used to extract such structure 
and parameters of the elastic scattering amplitude from the experimental data 
requires some different assumptions. If we assume that the hard pomeron is 
present in soft data, it will lead to a cross section of the order of 150 mb (a similar 
conclusion, via a different argument, has been reached in [13]). The uncertainty in 
this number is quite large, as unitarisation and saturation schemes are numerous. 
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Figure 7: The ratio of the real to the imaginary part of the amplitude (left), and 
the slope of the elastic differential cross section (right) as functions of t for the 
eikonalised amplitude at various energies: 100 GeV (plain curve), 500 GeV (long 
dashes), 5 TeV (short dashes) and 14 TeV (dash-dotted curve). 

Hence it seems that the total cross section can be anywhere between 108 mb [T6] 
and 150 mb. 

Other observables may be used to decide whether one has a simple extrapo- 
lation of the lower-energy data, such as in [TB] , or one is entering a new regime 
of unitarisation. Indeed, in the presence of the hard Pomeron, the saturation 
effects, which must then be present at LHC energies, can change the behavior 
of the real part of the cross section making it smaller than expected, especially 
in the near-forward region, and of the slope of the differential elastic scattering 
cross section. 

Despite the lack of an absolute prediction for total cross sections, the observa- 
tion of such features would be a clear sign that a new regime of strong interactions 
has been reached. 
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